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TECHNICAL MEMORANDUM 


INTERIM REPORT ON Ml CROFISSURING OF INCONEL 718 

INTRODUCTION 


When Inconel 718, a major constituent of the Space Shuttle Main Engine, is 
welded, microfissures are often found between the grains of the inetal in the heat- 
affected zone (HAZ) of the weld. Microfissures were discovered in Inconel 718 in 
1965 [1]. Since that time a great amount of effort devoted to studying the metallo- 
graphical setting in which microfissuring occurs failed to yield an understanding 
adequate to control it, so attention was turned toward the dynamics of the micro- 
fissuring process itself [ 2] . 

In 1976 [3], Professor Y. A rata and t' team of Japanese researchers from 
Kawasaki Heavy Industries, Ltd., publishea microfissuring envelopes separating 
regions of higher weld speed /electron beam (BB) power settings at which micro- 
fissuring was encountered in EB welds from lower speed /power settings at which 
microfissuring did not occur. They covered speeds to 60 inches per minute (IPM) 
and powers to 7 kW. Microfissuring envelopes for Hastelloy X in several heat treat 
conditions, as well as for Incoloy 800 and Incoloy 901, were plotted. 

The microfissuring envelopes of Arata et al. comprise a potentially very useful 
tool not only for eliminating microfissuring through adjustment of processing parame- 
ters but for characterizing the microfissuring tendency of a particular metal as well. 

Microfissuring appears to be associated with grain boundary liquation [4,5] 
below the solidus temperature, but microfissuring does not necessarily occur when 
grain boundary liquation is present [6]. This implies that grain boundary liquation 
can set up a condition of brittle vulnerability at temperatures slightly below the 
alloy solidus, but a second condition, sufficient mechanical strain, must also be met 
before microfissures form. Cracks tend to open up in a direction perpendicular to the 
the maximum tensile strain. A study of the orientation of microfissures [7] suggests 
that the strains responsible for microffisuring are oriented in a direction perpendi- 
cular to the plane of the workpiece surface or (often the same thing) parallel to 
the weld puddle boundary. 

Large grains, produced at temperatures above 1825°F at which the grain- 
boundary-pinning delta phase dissolves, appear to be a necessary but not sufficient 
condition for microfissuring [8,9]. The effect of large grains could be to concentrate 
detrimental impurities in reduced grain boundary area. A simple mechanical explana- 
tion for the grain boundary effect on microfissuring, however, comprises a limited 
deformation capacity for a single grain boundary so that a strain concentrated at 
liquated grain boundaries becomes more potent the fewer the grain boundaries avail- 
able to absorb it. The latter explanation is more attractive because of its simplicity 
and because it is easily quantified. As will be seen further in this report, the 
latter explanation, when quantified, appears to fit the data available at present. 


It is possible to build not merely a plausible qualitative picture, but a rough 
quantitative model of the microfissuring process from the hbove clues and a physical 
understanding of .the weld process. This has been accomplished and, pending further 
experimental study, the model appears to fit the experimental data and to offer 
predictive capability. 


TEST PROCEDURE 


Eight microfissuring specimens of 0.25-in. thick Inconel 718 plate were cut 4-in. 
wide by 12-in. long. Four of the plates were heat treated for 1 hr at 2150°F followed 
by an air cool. 

The object in heat treating half the specimens was to make them more prone to 
microfissuring. Large grains appear to be associated with microfissuring. Figure 1 
shows the results of grain size measurements taken on Inconel 718 samples exposed to 
1750, 1850, 1950, and 2150°F for 1 hr followed by an air cool. An ASTM grain size 
of 6 or 7 is normal for as received wrought Inconel 718 plate. An ASTM grain size 
around 2 would be expected to result from the above heat treatment. 
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Figure 1. Grain growth in Inconel 718 during 1 hr exposure to various 
temperatures followed by air cool to room temperature. 
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An EB bead-on-plate weld was run down the center of each plate using the 
voltage /current /speed parameters listed in Table 1. As shown in Figure 2 closely 
equivalent weld profiles were obtained with some focus adjustment. A mean weld root 
of 0.099 in. with a standard deviation of 0.010 in. and a mean weld crown of 0.137 in. 
with a standard deviation of 0.005 in. were obtained. A tendency of the weld fusion 
metal to drop down below the plate surface was observed, most pronounced at the 
slowest speeds (specimens 4 and 8). 


TABLE 1. WELDING PARAMETERS USED 


Specimen No. tt 

Voltage (kV) 

Current (mA) 

Speed (IPM) b 

1.5 

140 

40 

40 

2,6 

140 

30 

30 

3.7 

140 

20 

20 

4.8 

140 

10 

10 


a. Specimens 1 through 4 were heat treated at 2150°F for 1 hr 
followed by an air cool. Specimens 5 throuth 8 were as 
received. 


b. In all cases the energy delivered was 8.4 kJ/in. An increase in 
current (power) was compensated by an increase in weld speed. 



Figure 2. Crown and root widths of EB welds in 0.25-in. thick Inconel 718 
plate made with the same energy per unit weld length but at different 
weld speed/power conbinations. Focus adjustments needed to 
maintain approximately the same crown and root 
width are also shown. 
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The weld top surface was ground down O.OjO in. below the plate surface to 
reach the microflssured regions in steps of 0.010 in. The surface was etched [50 
percent HC1-50 percent HgOg (30 percent by volume)] after grinding. A fluorescent 

dye penetrant was then used to reveal the microfissures. 


TEST RESULTS 


Three of the four heat-treated specimens exhibited microfissures 0.030 in. below 
the top surface of the plate. The number of microfissures appeared to increase by 
roughly 2.25 microfissures per inch per kilowatt as shown in Figure 3. 

No microfissures were exhibited by the four as-received plates, which were 
subjected to the same weld procedures as the heat-treated plates. 


THEORY OF MICROFISSURING 


Microfissuring is thought to take place due to the action of tensile strains in the 
weld HAZ during the period when the HAZ begins to cool down and contract while its 
environment continues to heat up and expand. The tensile strains result as the 
shrinking HAZ is forced to conform to its still expanding environment. 

The mathematical model of microfissuring constructed here comprises three parts: 

1) A model of the temperature field in the vicinity of the weld. 

2) A computation of the stress-strain state in the vicinity of the weld due to the 
thermal stresses produced by the temperature field. 

3) A criterion for formation of microfissures under the stress-strain conditions 
prevailing in the weld HAZ. 

The model presented here employs only elementary components necessary to 
make an approximate calculation of the likelihood of observing microfissuring under a 
given weld situation. Experimental work is continuing to fully validate this model or 
an improved version. Improvements in the model are feasible and may turn out to be 
necessary in the light of future experimental data. 

The temperature field around the weld is modeled by an adaptation of the 
Rosenthal moving heat source model [10]. A moving point source on the surface of a 
plate of finite, specified width is combined with a moving line source extending down 
from the point source through the plate. The line source represents the spike of an 
EB weld; the point source, the nailhead, i.e., the extended region at the top of the 
weld, conceived as due to lower power density edge portions of the electron beam 
which fail to penetrate deeply and are absorbed at the top of the weld. 

The temperature rise T-T q above ambient temperature T q due to a point source 

of power P. moving over the surface of the plate in the X -direction with velocity V 
is given by: 
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gure 3. Microfissures per inch versus beam power and approximate 
location of micro fissu ring envelope point on speed /power 
constant weld cross section line. 
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T-T = 


p i 

T3< 


J5 < rtx > 


♦z 

n=l 
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h <V X> 
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V 

25 


< r n’«) 


n 


n 
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where 


k = thermal conductivity of metal 
a = thermal diffusivity of metal 

r n = /x^ + y^ + (z - 2nd)^ 

r ~ 9 o o 

r n * = / x* + y + (z + 2nd) z 


( 1 ) 


x,y,z = location of temperature site with respect to point source at 0,0,0. x 
is in the direction of motion; y is directed to the side; and z, is 
directed down into the interior of the plate perpendicular to the plate 
surface 

d = plate thickness. 


The summation, which is, of course, truncated in the computer calculation, 
comprises image sources which produce the effect of the back surface of the plate. 

The temperature rise T-T q above ambient due to a moving line source of total 
power P 2 distributed over the plate thickness d is given by: 


T-T 





where 

R= /77y T 

K q = zero order modified Bessel function of the second kind. 


( 2 ) 


An EB weld can be characterized by specifying Pj and In this work two 

"loss" parameters are used: Fj, the overall loss parameter, given as a percent and 

defined such that 


P 2 



( 3 ) 
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where P = total beam power; and F 0 . the nailhead loss parameter, given as a percent 

O a 

and defined such that 


= rdo 1 (P i + p 2 > 
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( 4 ) 


These losses can be calculated by matching the crown and root widths of a weld 
with the computed maximum widths between the solidus temperature isotherms at the 
top and bottom of the weld. The results of such calculations for the experimental 
welds of this report are shown in Figure 4. 

The difficulties in calculating the stress-strain history in the weld HAZ were 
circumvented by a crude but highly simplifying approximation. The strain imposed 
upon the weld HAZ is assumed equal to the net difference in thermal expansion 



Figure 4. Overall and nailhead power losses versus weld speed for constant 
energy per length equivalent profile KB welds. 
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between the HAZ itself and its "environment. " The "environment" is tentatively 
assigned the temperature encountered at a distance twice as far out from the weld 
centerline as the width of the weld. Thus given thermal expansion coefficient y . 
strain increment A £ is encountered such that 


Ae = p (AT ft - AT haz ) (5) 


where 


AT e = environmental temperature increment 

* ATjiaz = HAZ temperature increment, as one moves from initially unstrained 

material back along the weld in the HAZ. 


If no plastic strain were to occur the strain increment would be accompanied by 
a stress increment Ac, where 


Ao = E Ac 


( 6 ) 


where E = the linear elastic modulus of the metal. The use of a simple linear stress- 
strain relation is believed justified partly by the presence of molten or very soft 
material in the fusion zone of the weld preventing triaxial stress conditions during the 
imposition of a significant part of the thermal straining and partly by the inherent 
level of approximation which would make greater precision of one isolated aspect of the 
stress calculation superfluous. 

A nominal stress is calculated by adding An to its initial value a. 


o v , = o + Ao 


(7) 


and compared with an estimated flow stress o. t which vanishes at the metal melting 
point , 


°f ' K s (T M ' T HAZ ) 


( 8 ) 


where 


K = constant 
s 

Tm = melting (solidus) temperature of metal 
Thaz = HAZ temperature. 
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If the nominal stress a N would go outside the elastic region ”°f < ° n < + o r f» 
then a plastic strain increment Ae^ is computed sufficient to relax the nominal strain 
to a value that can be supported by the material. Thus, if 

a N > °f * 

then 

A £p = A £ * F l 
and a new stress is c 

C = Oj. 

If a N < - Op t hen 

• 

and 

a = -Oj. (12) 


(V)] ■ 


’omputed 


( 9 ) 


( 10 ) 


Otherwise , 


A Sp =0 . 


( 13 ) 


and 


0 = 0 , 


N ‘ 


Lastly, a damage theory of mi crofissuring is invoked, i.e. , microfissuring is 
assumed to occur when 


. A e„ 


( 14 ) 
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where Ac p / e c is the ratio of plastic strain increment to the critical plastic strain that 

would cause fracture under conditions of fixed temperature. Each plastic strain incre- 
ment is assumed to contribute to the final fracture condition in proportion to its ratio 
to the critical fracture strain. This assumption must be regarded as tentative and 
subject to revision as fu* + her data is acquired. The precedent for use of a damage 
concept comes from fatigi research. As in the case of fatigue there is no reason 
to use a linear damage criterion except for simplicity's sake in the absence of experi- 
mental or theoretical justification for a more complex damage relation. 

A critical fracture stress relation of the form 


e = t: 

C CO 


(T-T q ) 2 

B *~ 

e 



(15) 


where e^, T G , B are constants, and GS = grain size of metal, is assumed. 

R. G. Thompson's plot [11,12] of temperature versus critical strain for a hot 
tensile test of Inconel 718 with an ASTM grain size of 4 supplied values for the 
constants : 


e = 0.002 
co 


T C = T M * 25 ° F 


B = 2058 (°F) ‘ 


The grain size dependence assumes that a given extension per grain boundary Ad G 
is required for fracture. If the grain diameter is d G then ’ 



(16) 


By definition the number of grains per square inch at a magnification of xlOO is 

„GS - 1 , 

2 , so that 



( 17 ) 
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' CO 

( "') (iS 4 


tiS/2 




( 18 ) 


This accounts for the grain size correct ion exhibited in equation (15). 

These mathematical pieces have lieen put together in a computer model written 
out Appendix 1. The model computes damage curves versus speed or power for 
selected loss conditions. A damage value in excess of one implies micro Assuring. 

Let it be noted in conclusion that a more detailed model should also take into 
account a hysteresis effect due to time variation in grain boundary wetting. This has 
not been incorporated into the present model, nor is it likely to be until much more 
experimental data on the miorofissuring process becomes available. 


COMPARISON OP THEORY AND EXPERIMENT 


With the assumptions incorporated into the theory of miorofissuring presented 
here, it was anticipated that some adjustment of tlx> model would be required to reach 
good agreement with empirical data. Nevertlieless, at this point , theory and experi 
ment appear to agree well. 

'fable 2 shows computer model predictions of power speed combinations at which 
miorofissuring begins. Table :l shows the excess of power over that which should, 
according to the computer model, produce mierofissuring. Where there is an excess 
of power, miorofissuring would be expected. Where there is an insufficiency of power, 
miorofissuring would not be expected. 

TABLE 2. COMP l IT KR MODEL PREDICTIONS OP POWER SPEED COMBINATION 

AT WHICH MICROEISSHRINC. BEGINS 


Specimen No. 

Estimated GS 

Power (kW) 

Speed (1PM) 

1 

•) 

1.7 

7 

•) 

*> 

1.7 

7 

a 

•) 

1.2 

6 

4 

*) 

1 . 1 

6 

5 

() 

a 

a 

(> 

6 

a 

a 

7 

r 

a 

a 

8 

(> 

a 

a 


a. 


Powers substantially above 20 kW and speeds substantially 
above 200 1PM would be required. 
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TABLE 3. EXCESS POWER ABOVE COMPUTER MODEL POWER 
REQUIREMENT TO INITIATE MICROFISSURING 


Specimen No. 

Excess Power 

1 

+3.0 

2 

+2.5 

3 

+1. 6 

4 

+0.3 

5 

a 

6 

a 

7 

a 

8 

a 


a. Substantially negative excess 
power (>-25 kW) predicts no 
micro Assuring. 


Figure 5 is a plot of measured microfissures per inch versus calculated micro- 
fissuring power excess for the four heat-treated specimens. Microfissuring only 
occurred for these specimens, where there was a calculated excess power. From the 
limited data it would appear that the amount of microfissuring increases with excess 
power at approximately 4.5 microfissures per inch per excess kilowatt. 



MICROFISSURING POWER EXCESS 
(KW - CALCULATED) 


Figure 5. Microfissures per inch versus microfissuring power excess. 
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CONCLUSIONS AND RECOMMENDATIONS 


Microfissures have been produced in 0.25-in. thick Inconel 718 plates by 
electron beam welding. 

A computer model of microfissuring has been worked out. The model is com- 
patible with the experimental data obtained so far. The power level required by the 
model to produce microfissuring can be computed for a given weld situation. About 
4.5 mierofissures per inch per excess kilowatt above the computed level to produce 
microfissuring were produced. Slowing the weld speed and reducing power below the 
level required to produce mierofissures should eliminate microfissuring. 

The model incorporates effects of weld power and speed, weld loss or efficiency 
parameters, and grain size. It can be adapted to other metals subject to micro- 
fissuring by altering material parameters. 

Further empirical work, particularly on thicker plate at higher powers, is 
needed to provide a more substantial data base suitable for thorough evaluation of the 
model . 


Once a means for producing microl'issured and unmicrofissured , but otherwise 
similar, welds has been perfected, a study of t lie effect of microfissuring on mechanical 
properties can and slnmld be carried out. Some mechanical testing using the speci- 
mens prepared in the previously eited study is now being planned. 

Further editing of the computer program to improve style and efficiency is 
needed. It is planned to do this in the course of preparing this and other weld 
model programs for public distribution through the Computer Software Management 
and Information Center at the University of Georgia. 
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APPENDIX A 

MICRO FISSURING MODEL COMPUTER PROGRAM LISTING 


Overall and nailhead losses plus the desired weld root width are inserted into 
the program. A velocity range, usually 0 to 100 inches per minute, is preselected 
within the program but can be changed if desired. 

The power settings corresponding to the velocities are computed as well as the 
location of the peak temperature position of the HAZ. Starting at the beginning of 
cool-off the plastic strain associated with the thermal stresses is estimated for a 
series of positions from the cool-off point to a point 2.5 in. beiand the weld. A 
nonlinear scale weighted to take more points at higher temperatures is used. 

From the difference between the estimated critical fracture strain (equation 15) 
and the plastic strain calculated as explained above is calculated a damage parameter, 
which is plotted versus power and velocity as program output. If the damage 
parameter exceeds unity, micrufissuring is predicted. 

The program is written in TEK SPS VO 2- 01 Basic for use on a Tektronix CP 
1164 (DEC PDP 1134) central processor incorporated into a Tektronix WP 1220 Mod 04 
Data Acquisition System. The pi*ogram runs in a little over an hour. 


10 F'kM ♦** MICPOFISSURIMG DAMAGE calculation *** 

.*« SOI »♦» DATA input *** 

'50 H LEAP 

4u PRINT "PLEASE SELECT THE NUMBER OF THE METAL TO BE MELDED" 

SO PRINT “FROM THE LIST BELOW "\PRINTsPRINT\PRINT 

60 PRINT "0 UNLISTED METAL "PRINT 

TO PRINT "1 2219 ALUMINUM" 'PRINT 

80 PRINT "2 302 STAINLESS STEEL"'PRINT 

•?0 PRINT "3 321 STAINLESS STEEL "nPRINT 

100 PRINT "4 INCONEL 718" PRINT 

110 INPUT 0 

120 IF 0=1 THEN 460 

13M IF 0=2 THEN 490 

140 IF 0=3 THEN 520 

150 IF 0=4 THEN 550 

160 PRINT "WHAT IS THE METAL TO BE WELDED? "vPR INT SPRINT \PRINT 
170 INPUT M» 

180 PRINT "PLEASE WRITE IN THE LOWER AND UPPER LIMITS OF THE MELTING" 
190 PRINT "TEMPERATURE RANGE OF THE ALLOY TO BE WELDED USE UNITS OF" 
20O PRINT "DEGREES FAHRENHEIT SEPARATE THE TWO UALUES BY A COMMA." 

210 PRINT "NOTE A G000 SOURCE REFERENCE FOR THIS KIND OF DATA IS THE" 
220 PRINT "AEROSPACE STRUCTURAL MATERIALS HANDBOOK "\PRIHTvPRINT\PRINT 
2 30 INPUT LT.IJT 

240 PRINT "WHAT IS THE THERMAL CONDUCTIUITY OF THE ALLOY IN BTU’S PER" 
250 PRINT "HOUR PEP FOOT PER DEGREE FAHRENHEIT?"\PRINTnPRINTnPRINT 
260 INPUT K1 
270 LET F1=KU2 441E-05 

280 PRINT "WHAT IS THE SPECIFIC HEAT OF THE ALLOY IN BTU’S PER LB" 

290 PRINT "PER OEGPEE FAHRENHEIT?"\PRINTsPRINTnPRINT 
300 INPUT Cl 

310 PRINT "WHAT IS THE DENSITY OF THE ALLOY IN LBS PER CUBIC INCH?" 

320 PRINT PRINT PRINT 
330 INPUT RO 

340 LET A1*K1^R0/CU56 83 
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350 HrtfkT "WHAT IS THE COEFFICIENT OF THERMAL EXPANSION • 

360 PRINT “IN INCHES PER INCH PER DEGREE FAHRENHEIT?" 

370 Iff’UT MU 

380 PRINT “WHAT IS THE ELASTIC NOOULUS IN PSI?" 

390 INPUT E 

400 PRINT "ASSUMING THAT. IN THE NEIGHBORHOOD OF MELTING. TIC" 

410 PRINT “FLOW STRESS IS PROPORTIONAL TO THE SQUARE OF THE " 

420 PRINT -DIFFERENCE BETWEEN HAZ TEMPERATURE AND MELTING " 

430 PRINT -TEMPERATURE. WHAT IS THE FLOW STRESS COEFFICIENT" 

440 PRINT "IN PSI PER DEGREES FAHRENHEIT SQUARED?" 

450 INPUT KS 

460 LET UT-1190VLET LT-lOIOsLET Kl-2 44E-83VLET A 1-5 56 

470 LET MU-1 44E-05 -LET E-SOOOOOevLET KS* 08 

480 LET MS--2219 ALUMINUM" ''GOTO 588 

490 LET UT-2650NLET LT-2550VLET Ki-4.Z7E-04vLET AI-.5 

500 LET MU-1.2E-05SLET E«* 6E+07MXT KS- 05 

510 LET MF--302 STAINLESS STEEL” ''GOTO 588 

520 LET UT-2558MET LT-2580M£T KI-4.83E-84MET Al- 66 

530 LET HU-1 15E-85MET E-l 6E*07 nLET KS- 85 

540 LET NS--321 STAINLESS STEEL -'GOTO 588 

550 LET UT-2437''LET LT»23O0\LET K1-3.83E-84MET Al- 71 

560 LET MU-9 8E-06M.ET E-l 5E407MET KS- 856 

570 LET N$-" INCONEL 718"nG0T0 588 

580 PR 1NTPR I NT\P* I NT 

590 PRINT "WHAT IS THE ASTM GRAIN SIZE OF THE METAL?" 

680 PRINT vPRINTnPRINT 
610 INPUT GS 

620 PRINT "THE METAL TO BE WELDED IS NON CHARACTERIZED. ME NON" 

630 PRINT "TURN TO THE WELDING PROCESS PARAMETERS . "nPRINTnPRINTnPRINT 
640 PRINT "WHAT IS THE AMBIENT TEMPERATURE OF THE METAL IN DEGREES" 

650 PRINT "FAHRENHEIT? IF THE METAL IS PREHEATED. GIUE THE PREHEAT" 

660 PRINT "TEMPERATURE AS THE AMBIENT TEMPERATURE. "nPRINTnPRINTnPRINT 
670 INPUT TO 

680 PRINT "PLEASE SPECIFY THE PLATE THICKNESS IN INCHES." 

690 PR I NT' >PR I NT \PR I NT 
700 INPUT Ml 

710 PRINT "PLEASE SPECIFY WELD ROOT WIDTH IN INCHES " 

72C PR INT'PRINT 'PRINT 
730 INPUT NR 

740 PRINT "WHAT PERCENT OF TOTAL BEAN POWER IS LOST FROM THE HELD" 

750 PRINT "PUO0LE DUE TO PROCESSES OTHER THAN CONDUCT I OH BY THE PLATE?" 
760 PRINT "THESE LOSSES INCLUOE RADIATION ANO METAL EVAPORATION FRON" 
770 PRINT "THE VICINITY OF THE PUOOLE "nPRINTnPRINTnPRINT 
780 INPUT FI 

790 PRINT “WHAT PERCENT OF THE REMAINING BEAM POWER IS ABSORBED CLOSE" 
800 PPINT “TO THE METAL SURFACE SO AS TO FORM THE E8 HELD NAILHEAD?" 

810 PRINT “THIS WOULD BE THE PERCENTAGE OF THE BEAM CURRENT LACKING" 

820 PPINT "SUFFICIENT POWER DENSITY TO VAPORIZE THE METAL " 

830 FRINT' PRINT' PRINT 
840 INPUT F2 

850 PRINT "PLEASE SELECT POWER RANGE UPPER LIMIT IN KILOWATTS' " 

860 PPINT "5. 10. 25. OR 58." 

670 INPUT KW 

880 PRINT "TO PROCEED PRESS THE SPACE BAR " 

890 WAIT 
900 PAGE 
910 ONEPP NOMARN 

920 REM *** CONSTANT PENETRATION LINE PLOT *** 

930 DELETE "A.Ph 
940 LET M-20 

950 OIM UA< M > • PA< M > » XS( M ) 

960 LET "7*100 

970 FOR 1*8 TO M 

988 LET UA< I >*UZ»' < !♦! V* M*1 ) ) 

990 LET "1-UA' I > 

1808 G06U6 2340 
101O LET XSC I >*X 
1020 LET PA< I >«P0 
1030 »CXT I 
1040 DELETE PB.'JB 
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BEAM POWER (KILOWATTS)" 


050 OIM PB(M>*Ue<M> 

060 FOR 1*0 TO M 
070 LET U6< I >-UA< I ) 

080 LET P8< I >*PA< 1 1 
090 IF PAtlXKM THEN 1130 

100 IF PA< I »*KM THEN IF PA< I-l »-KH THEN 1120 
110 LET U8< I >«UA< I-l VKUZAM^l ))*<<KW-PAC 1-1 >V(PAC !>-«* I-l >>> 

120 IF Pf* I »*KM THEN LET P0( I >-KH 
130 NEXT I 
140 PAGE 
150 GOSUB 1690 

160 FOR 1*0 TO IOnPRINT'-NEXT I 
170 PRINT "HELD SPEED** 

180 PRINT - ( IPM >** 

190 FOR 1-0 TO 14vPRINT\NEXT I 
200 PRINT - 
210 HIHOOH 0,KW,0,UZ 
220 UIEUPORT 260,763,156,650 

230 SETGR UIEM. WINDOW, TICS 10, 10,5,5.GRAT 6, 6,3, 3 
240 XYPLOT P8,«J0 
250 PRINT --C-W-C-L- 
260 DELETE DA 
270 DIM Qm(M> 

280 FOR J-0 TO M 
290 LET P2*PA( J»< l-F1^100) 

300 LET P1*P2*(F2'100> 

310 LET P2-P2-P1 
320 LET U1»UA<J> 

330 GOSUB 1760 
340 LET 0A(J>«0M 
350 NEXT J 
360 OELETE 06 
370 DIM OB<M> 

380 FOR 1*0 TO M 
390 LET 00< I >-0A< I > 

400 IF 0*. I »*2 THEN LET 0B( I >*2 
410 NEXT I 
420 PAGE 
430 GOSUB 1698 

440 FOR 1*0 TO 10\PRINT\NEXT I 
450 PRINT "HELD SPEEO" 

460 PRINT " < IPM >** 

470 FOR 1=0 TO 14\PRINTsNEXT I 
480 PRINT ** DAMAGE ** 

490 WINDOW 0.2,O,U2 
500 UIEUPORT 260,763,156,650 

510 SETGR UIEH.MINOOM.TICS 10, 10,5,5.GRAT 6,6,3,3 
520 XYPLOT 00, U6 
530 PRINT •••C'M^C^L** 

540 GOSUB 1690 

550 FOR 1*0 TO 10nPRINT\NEXT I 
560 PRINT "DAMAGE" 

570 FOR 1=0 TO ISnPRINWCXT I 
580 PRINT " 

590 WINDOW G. KM, 0,2 
600 UIEUPORT 260,763,156,650 

610 SETGR UIEM,MINOOH* TICS 10, 10,5,5, GRAT 6,6, 3, 3 
620 XYPLOT PBOB 

630 REM Nt DELAY TO PREUENT SKIPPING PLOT *** 

640 FOR 1*0 TO 1000 
650 LET B-EXP* 10 > 

660 ICXT I 

670 PRINT "•I'M 'L M 
680 ENO 

690 REM «»* CUR' € LABEL SUBROUTINE *** 

700 PPINT M*, **--", 

710 PRINT "PLATE THICKNESS ", HI." INCHES" 

720 PRINT **HA2 ";Y»" INCHES OFF CTR. ",Z;“ INCHES DEPTH G8- "iGS 
730 PRINT FI.-'. NONCONDUCT IUE LOSSES — " ■FZ.—'i TO NAILHEAO" 

740 PPINT MR," INCH ROOT UIOTH" 
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1750 RETURN 

1760 REM 111 DAMAGE CALCULATION SUBROUTINE *** 

1770 REM 111 COMPUTATION OP HA2 TEMPERATURE HISTORY *** 

1780 LET Z>.03sREM DEPTH OF MICROFISSURINC IS *98 INCHES 
1790 Q OSU B 2940 

1000 DELETE XA/XB/TA»TB/TC/TO/S»£F#EP»ET 
1810 LET N>288 
1820 LET XF-2 5 

1830 DIM NA<N)/TA<N)/T8<N)/TC<N)/T0<N) 

1848 LET TA-UTnLET TB-LT 
1850 FOR 1-0 TO N 

I860 LET XACI^XSCJlXXF-XSOJlttMOttMO 
1870 LET X— XA< I ) 

1880 G0SU8 2780 
1890 LET TC< I >*M2 
1900 LET Yl-Y 
1910 LET Y»2»Y 
1920 GOSUB 2788 
1930 LET TtX I )-M2 
1940 LET Y-Yl 
1950 NEXT I 

I960 REM »» COMPUTATION OF HAZ PLASTIC STRAIN HISTORY *** 

1978 OIM S<N)/EP<N> 

I960 LET S-0M.ET EP-0 
1990 FOR 1*1 TO N 

2000 LET OE-MUK T0< I >-TQf. 1-1 >-TC( I HTCC 1-1 ) ) 

2010 LET OS-ESOE 

2020 IF SGN< OS X >SGN< S< I~l > > THEN LET EP< 1-1 )-0 

2030 IF <LT-TC(I>»0 THEN LET SF-KS*< LT-TCCI > >-2 

2040 IF <LT-TC<I >><-e THEN LET SF-0 

2050 LET SN-S< I-l >*0S 

2060 IF SN>SF THEN 2090 

2070 IF SN<-SF THEN 2180 

2080 IF SN<»SF THEN IF SN>»- SF THEN 2110 

2090 LET £P< I XP( 1-1 )*OEt<. 1-<SF-S< I-l > VOS >\LET S< I >-*SF\GOTO 2120 

2100 LET EPC I >*EP< 1-1 VOE*' l+< SF+S* 1-1 > VOS VLET S< I >— SFM30T0 2120 
2110 LET SCI >-SN 
2120 NEXT I 
2130 DIM ET<N> 

2140 FOP 1*0 TO N 

2150 IF EPv I »-0 THEN LET ET< I >*EP< I > 

2160 IF EPCIX0 THEN LET ET<I >-0 
2170 NEXT I 

2180 REM l« FRACTURE ENUELOPE REPRESENTATION »U 
2190 DIM EF< N * 

2200 LET TOLT-25 
2210 FOR 1*0 TO N 

2220 LET EF< I >»2E-03*< 2'< < GS-4 V2 ) ) 

2230 LET TF*( TC< I >-TG >♦«. TC< I >-TG V2056 
2240 IF TF >88 THEN LET TF-83 
2250 LET EF< I >*EF< I >*EXP< TF > 

2260 NEXT I 

22 70 REM »« DAMAGE CALCULATION Ml 
2288 LET 0M*8 
2290 FOR 1-1 TO N 

2300 IF EF< I »9 99999E+34 THEN 2320 

2310 IF ET< I >>ET< 1-1 ) THEN OM-OM*ET< I VETU-1 >V«£F< I >»€F< I-l ) V2> 

2328 NEXT I 
2336 RETURN 

2340 REM ttt SUBROUTINE TO COMPUTE BEAM POWER *** 

2358 LET Y«HR^2 

2368 LET Z-Mi » 

2378 LET Tl-LT 

2380 LET X*0 

2390 LET OX— l 

2400 LET P0*1 

2410 LET P1-P0*<1-F1/100> 

2420 LET P2"PI*< 1-F2/188) 

2430 LET P1*P1-P2 
2440 GOSUB 2788 
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2450 LET M4-M2 ^ ^OOR 

2460 LET X-X+OX 
2470 GOSUB 2780 
2488 LET M5-M2 
2498 IF H5>M4 THEN GOTO 2498 
2508 IF M5-M4 THEN GOTO 2520 
2510 LET DX*< -OX >'10\GOTO 2450 
2520 LET P0-P0*< < T1-T0 V( N2-T0 ) ) 

2530 RETURN 

2540 REN m LOCATION OF HEAT -AFFECTED-ZONE *$* 

2550 LET Y* . l 'GOSUB 2670 
2560 LET OY*. I 
2570 LET N6>M2 
2580 LET Y-Y+OY 
2598 G08UB 2670 

2600 IF ABS<M2-LTX.ei*LT THEN 2660 
2610 IF H2<LT THEN IF N6<H2 THEN 2570 

2620 IF N2>LT THEN IF H6>N2 THEN 2578 

2630 IF M2<LT THEN IF H6>N2 T»CN LET OYX-OYVIO 

2640 IF H2>LT THEN IF N6<N2 THEN LET OY*<-OY>/10 

2650 GOTO 2S70 
2660 RETURN 

2670 REM m TEMPERATURE MAXIMI2AT10N SUBROUTINE *** 

2680 LET X-0 

2690 LET OX— .1 

2700 GOSUB 2700 

2710 LET N4*N2 

2720 LET X-X*OK 

2730 GOSU0 2780 

2740 IF M4<M2 THEN 2710 

2750 IF M4-M2 THEN 2770 

2760 LET OX— OX'IO'-GOTO 2710 

2770 RETURN 

2780 REM *** SUBROUTINE TO COMPUTE TEMPERATURE OF HEAT SOURCE ARRAY *** 
2790 LET Z2*Z 

2800 LET M2*0 
2810 LET 03*8 
2328 LET D 2*2 *141 
2830 GOSUB 3170 
2840 GOSUB 3000 
2850 LET M2*M1+T0*T2 
2860 LET M3=N2 
2870 LET D3-03+02 
2888 GOSUB 2920 

2890 IF < M2-M3 X 01*M2 THEN GOTO 2910 
2900 GOTO 2860 
2910 RETURN 

292© REM SUBROUTINE TO POO NEXT TWO HEAT SOURCES *** 

2930 LET 22*03*2 
294© GOSUB 3170 
2950 LET M2*M2N11 
2960 LET 22*03-2 
2970 GOSUB 3170 
2980 LET M2*M2N1l 
2990 RETURN 

30O0 REM *** SUBROUTINE TO COMPUTE LINE HEAT SOURCE TEMPERATURE 
3O10 REM DISTRIBUTION *** 

3028 LET K0-0 

3038 LET T2*0 

3040 LET AP*U1-'2'A1»X 

3050 LET Z3*Ul/2/Ai*SGR<X*X*Y*Y> 

3060 IF AR/88 THEN 3160 

3070 GOSUB 3250 

3060 IF AR<0 THEN 3110 

3090 LET T2*P2'U1/6.2032''K1>'EXF<:aR>*K0 

3100 GOTO 3160 

3110 LET AR— AR 

3120 IF 23 >20 THEN 3150 

3130 LET T2-R2/HI/6 2832'K1*EXP< AR >*K0 

3140 GOTO 3160 


Rage & 
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3150 LET T2-P2/W1 '6 2832/RUl 2S33/EXP<Z3*«R>/90R(Z3> 

3160 RETURN 

31 70 REN W SUBROUTINE TO COMPUTE TEMPERATURE tU 
3160 LET Sl-SQR<X*Xm , i>Z2*Z2> 

3190 IF $1-0 THEM LET Ml-UUIBOOMF Sl-0 THEN 3240 
3200 LET AR-U1- 2/Al*(Sl«tt> 

3210 IF hR> 88 THEN LET Ml^ 

3220 IF NR >88 THEN GOTO 3240 

3230 LET M1*4>1 '2'3 1416'KI/S1/EXP<AR> 

3240 RETURN 

3250 REM *#* HOOIFIED BESSEL FUNCTION, SECONO Ki 0. ZEROTH ORDER lit 
3260 REM II* POLYNOMIAL APPROXIMATION *#* 

3270 IF Z3>88 THEN 3400 
3280 IF Z3>2 THEN 3350 
3290 GOSUB 3420 
3300 LET G2-Z3SZ3/4 

3310 LET K0-* < 7 . 4E-06FG2+1 . 073E-04 >tG2*2 . 62698E-03 >HG2 

3320 LET KO-C C K0» 0348859 >*G2^ 2306976 >*G2 

3330 LET K0X<KO* 4227842 >*G2- 5772157)- 5*LOG<G2>*I0 

3340 GOTO 3410 

3350 LET G2-2/Z3 

3360 LET K0X(5 32O8E-04FG2-2 5154E-03 >*G2*5 87872E-03 >*G2 

3379 LET KOX<K0- 01062446 )*C2* 02189568 )*G2 

3388 LET KO-UKO- 07832358 XC2*1 253314 VSQR< Z3 )/EXP(Z3> 

3390 GOTO 3410 
3400 LET K0-0 
3410 RETURN 

3420 REM *** HOOIFIED BESSEL FUNCTION, FIRST KINO. ZEROTH OROER «** 
3430 REM *** POLYNOMIAL APPROXIMATION *** 

3440 IF Z3>3 75 THEN 3S00 
3450 LET Gl-Z3*Z3 

3460 LET I0«<5 923979E-10»G1*6 S6817E-08>*G1*6 80123E-06)*G1 
3470 LET I0-(<10*4 3394E-04>*G1 4- 8156252 >*C1 
3460 LET I0»< 104- 23 )*G1+1 
3498 GOTO 3540 

3500 LET I0-U153 445/23-171 822 VZ3*73 2919VZ3 

3510 LET 10** t 10-15 2595VZ3M 81198>'Z3 

3520 LET I0*<< 10- 0830909 >Z3* 0316855 VZ3 

3538 LET I0*< 1 10* 0496222 > Z3^ 3989423 )*EXP<Z3>/S0R<Z3> 

3540 RETUPN 

READY 

» 
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APPENDIX B 

SAMPLE MICROFISSURING MODEL OUTPUT 


The microflssuring model computer program prints: 

1) A curve of weld speed versus beam power for the assumed loss parameters. 

2) A curve of damage versus beam power. 

3) A curve of weld speed versus damage. 

The weld speeds and beam powers for which the damage exceeds unity are 
taken to exhibit microflssuring. 


INCONEL 7 ! 8 — PLATE THICKNESS 23 INCHES 
HMZ 836 INCHES OFF CTR 23 INCHES DEPTH GS- 2 
20 3\ HONCOHOUCTIUE LOSSES— 7 55; TO NAILHEAO 
1 12 INCH ROOT HIDTH 


HELD SPEED 
»' IPH> 



e 


2 4 « 


13 3 7 

SEAN POWER <KlLOUSrre> 




9 


It 


Figure 6. Computer output : weld speed versus beam power. 
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INCONEL 718 — PLATE THICKNESS 29 INCHES 
HHZ 069 INCHES OFF CTR 03 INCHES DEPTH. CS- 2 
20 NONCONDUCT IUE LOSSES— 7 33; TO NAILHEAO 

112 INCH ROOT WIDTH 


OANAGE 



0 2 4 6 8 10 

1 3 3 7 9 


BERN POWER (KILOWATTS) 

Figure 7. Computer output : damage versus beam power. 


INCONEL 718 — PLATE THICKNESS 25 INCHES 
mz 069 INCHES OFF CTR 03 INCHES DEPTH GS- 2 
2C 9?. NONCONOUCTIUE LOSSES— 7 3* TO NAILHERO 
112 INCH ROOT WIDTH 


HELD SPEED 
<IPH> 



0 


.4 


2 


.0 

6 1 
ORHRGE 


1.2 16 2 

1.4 1.0 


Figure 8. 


Computer output: 


weld speed versus damage. 



APPROVAL 


ORIGINAL PAGE IS 
OF POOR QUALITY 


INTERIM REPORT ON Ml CROFI SSURING OF INCONa718 

By A. C. Nunes, Jr. 


The information in this report has been reviewed for technical content. Review 
of any information concerning Department of Defense or nuclear energy activities or 
programs has been made by the MSFC Security Classification Officer, This report, 
in its entirety, has been determined to be unclassified. 



Director, Materials and 'Processes Laboratory 


eu.S. GOVERNMENT PRINTING OFFICE 19*3 1462-238 


23 


